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Abstract 

The Feynman diagram generator FeynArts and the computer algebra program 
FormCalc allow for an automatic computation of 2 — ► 2 and 2^3 scattering 
processes in High Energy Physics. We have extended this package by four new kine- 
matical routines and adapted one existing routine in order to accomodate also two- 
and three-body decays of massive particles. This makes it possible to compute auto- 
matically two- and three-body particle decay widths and decay energy distributions 
as well as resonant particle production within the Standard Model and the Minimal 
Supersymmetric Standard Model at the tree- and loop-level. The use of the program 
is illustrated with three standard examples: h — > bb, [i — ► eu e Va, and Z — ► v e v e . 
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1 Introduction 

The analytical calculation of Feynman diagrams and their subsequent numer- 
ical evaluation represent the standard technique for calculating cross sections 
of scattering processes and widths of particle decays in High Energy Physics. 
Unfortunately, the number of diagrams to be computed grows rapidly with 
the number of external particle legs, internal particle loops, and exchanged 
particles. This makes an automatic generation and computation of Feynman 
diagrams a highly desirable feature, not only to reduce the time needed for 
the calculations, but also to eliminate possible sources of errors. 

For more than ten years, a well-tested and easy-to-use Feynman diagram gen- 
erator has existed in the form of the Mathematica package FeynArts [1] . More 
recently, a fast and equally reliable computer algebra program for the eval- 
uation of the Dirac and color traces inherent in the generated diagrams has 
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been implemented within Form (FormCalc) [2]. Together with LoopTools [2] 
(an adapted version of the Fortran program FF [3] for the numerical evalua- 
tion of tensor loop integrals), the FeynArts and FormCalc program package 
allows for an automatic computation of tree- and loop-level 2 — * 2 and 2^3 
scattering processes within the Standard Model and the Minimal Supersym- 
metric Standard Model. Since the phase space for each particular scattering 
process depends on the number of in- and out-going external particles, it is 
best evaluated directly within Fortran with specific kinematical routines. So 
far, only 2^2 and 2^3 scattering processes have been implemented within 
FeynArts and FormCalc. The aim of this Paper to extend their applicability 
also to two- and three-body decays of massive particles and to resonant par- 
ticle production. We provide four new kinematical routines and one adapted 
routine, which allow for the computation of the total two- and three-body de- 
cay widths and of the energy distribution of one decay product for three-body 
decays. 

The remainder of this Paper is organized as follows: In Sec. 2 we collect the 
well-known kinematical formula? for two- and three-body decays. Their imple- 
mentation in the new and modified Fortran routines is described in Sec. 3. As 
a test of the new routines, we calculate in Sec. 4 three standard examples for 
scalar, fermion, and vector boson decays. A short summary is given in Sec. 5. 



2 Kinematics for two- and three-body decays 



The decay width of a heavy particle with mass m 1? 
dTi = j\M\ 2 d^i, 



depends on the flux factor F, i.e. the number of decaying particles per unit 
volume, the squared invariant scattering amplitude \A4\ 2 , and the phase space 
volume d$j for i final-state particles. In the rest frame of the decaying particle, 
k\ = (mi, 0, 0, 0), the flux factor is simply 

F = 2mi. (2) 



The four-momenta of two decay products with masses m 2 , m 3 , 
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are uniquely defined by their energies Ef = \ki\ +mq, squared three-momenta 
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and by the center-of-mass scattering angle 9. The phase space volume for 
two-body decays is then [4] 
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where the dependence on the azimuthal angle integrates to 2tc. The total 
two-body decay width is obtained by integrating over the polar angle 9. [Note 
that in the rest frame of the decaying particle, the distribution in the polar 
angle 9 depends on the arbitrary definition of a reference axis.] 

For three-particle decays, the phase space volume is [4] 



d$ 3 = , ^ , ^ , ^ (27r)^W(fc 1 - k 2 - k 3 - h) 
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Here, the 5-function has been used to eliminate the integration over k% and #2, 
and the trivial integration over the azimuthal angle 04 has already been per- 
formed. For three-body decays, we will not only consider total decay widths, 
but also distributions in the decay energy £4. 

While the standard version of FeynArts and FormCalc diverges for 2 — > 2 and 
2^3 cross sections, which proceed through an intermediate s-channel reso- 
nance, our modified version makes it possible to calculate also these resonance 
cross sections using the formula 

47rT(R^ab)T(R^X) 
a am (^s) - (s _ m|)2 + m|(r , t)2 . (7) 



where s is the squared center-of-mass energy of the initial state particles a and 
b, rriR is the mass of the resonance R, and the decay widths are calculated as 
described above. 
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3 Description of the new and modified Fortran routines 

The phase space integrals described in the previous Section have been imple- 
mented in FeynArts and FormCalc with two new Fortran routines for two- 
and three-body decays, 

• lto2.F, 

• lto3.F, 

and two corresponding common blocks, 

• lto2.h, 

• lto3.h. 

The new Fortran routines have been adapted from the existing routines for 
2^2 and 2^3 particle scattering by modifying variable and center-of-mass 
energy definitions, the flux factor, and the initial-state four-momenta. The 
Fortran routine 

• num.F, 

which handles the averaging and the summation over initial and final state 
spins, has been modified to accomodate single-particle initial states. For this 
purpose, the routine num . F checks if the new preprocessor variable DECAY has 
been defined in lto2.F or lto3.F or not.Q 



4 Three examples: h — > bb, jj, — > ei/ e z^, and Z — > v e v e 

In order to demonstrate that the new kinematical routines work properly 
within the extended package of FeynArts and FormCalc and are straightfor- 
ward to use, we compute in this Section three standard decay widths. Note 
that the FormCalc routine process .h has to be adapted to every considered 
scattering process. In particular, the correct phase space generator has to be 
included at the end of process .h. However, the calling structure of the main 
program run . F has not been changed. For more details we refer the reader to 
the FeynArts [1] and FormCalc [2] manuals. 



The new Fortran routines can be obtained from klasen@mail.desy.de. 



4 



777 u /OpV 


r(ff^hh) /MpV 


1UU 


A 1 K 
4.10 


i in 

11U 


4.5o 


1 on 


k nn 


130 


5.42 


140 


5.85 


150 


6.27 



Table 1 

Decay widths of a Standard Model Higgs boson into a pair of bottom quarks for six 
different Higgs masses ran- 



4.1 h^bb 



The isotropic, scalar two-body decay h — > bb represents the dominant Higgs 
decay mode for a Standard Model Higgs boson with mass below m H = 140 
GeV. Using the default parameter settings of FeynArts and FormCalc and the 
Standard Model initialization file sm_ini .F, we obtain the decay widths listed 
in Table 1. Identical results are obtained with the tree- level formula [5] 



r(H -> bb) 
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where 1/a = 137 is the inverse of the electromagnetic fine structure constant, 
sin 2 6\v = 1 — m^ v /m 2 z is the electroweak mixing angle, Nc = 3 accounts for 
the color degree of freedom, m#, mz = 91.1882 GeV, and mw = 80.419 GeV 
are the masses of the Higgs, Z-, and lU-bosons, and m b = 4.7 GeV is the mass 
of the bottom quark. We have checked that the cross section differential in 
the scattering angle 9 is indeed constant. 



4.2 /i -> ez/ e z/ M 



In the Standard Model, the three-body decay of muons with mass m M = 
105.658 MeV into electrons and neutrinos proceeds through a virtual W-boson 
and accounts for almost 100% of the muonic decay width. With the default 
parameter settings of FeynArts and FormCalc, we obtain 



T(fi -> ez/ e z/ M ) = 2.8 x 10~ 10 eV 



(9) 
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which agrees with the result obtained from the lowest order formula [6] 

G 2 m 5 

r> - eu e u,) = ^-f , (10) 



where we have averaged over the initial-state spin, Fermi's constant G has been 
computed from G = n a/(y/2 sin 6^ m 2 ^), and the electron mass has been 
neglected. Furthermore, the calculated muon life time r = h/Y = 2.3 x 10~ 6 s 
is in good agreement with the measured value of 2.2 x s. As mentioned 
above, for three-body decays we consider also the energy distribution of one 
of the decay products, in this case of the observed electron [6] 

dr , x G 2 ml _ / 4E e \ , N 



As can be seen in Fig. 1, the numerical evaluation of Eq. (11) agrees very well 
with the result of FeynArts and FormCalc. 



4.3 Z -> u e u e 



In our final example, we consider the invisible branching ratio of the Z-boson 
in the Standard Model, which represents an important boundary condition for 
the number of neutrinos realized in nature. The tree-level formula [6] 

Y(Z -> v e v e ) = 9 amz (12) 

V e! 24 sin 2 9 W cos 2 9 W V ; 



and our modified program package FeynArts and FormCalc unanimously lead 
to the numerical result 

Y(Z -> v e v e ) = 160 MeV. (13) 



For three generations of neutrinos and employing the measured total width of 
the Z-boson, r tot = 2.49 GeV, this translates into an invisible branching ratio 
of 

BR inv = 3 r(Z T T^ e) = 19.3%, (14) 

1 tot 



which is close to the experimental value of 20.0%. 



6 



x 10 
0.12 

0.1 

0.08 



^0.06 



0.04 



0.02 



-16 







Electron Energy Spectrum in \x — > ev e 



Eq. (11) 

FeynArts & FormCalc 




0.01 0.02 0.03 0.04 0.05 0.06 

E e [GeV] 



Fig. 1. Electron energy spectrum of the decay width for muons decaying into elec- 
trons and neutrinos. 



5 Summary 



In this Paper, we have described the implementation of two- and three-body 
decays and resonance cross sections in the computer algebra package FeynArts 
and FormCalc. We have provided four new Fortran routines, which generate 
the phase space for two- and three-particle decays, and we have modified one 
existing routine in order to account for correct spin averages of one-particle 
initial states. While the new routines have been obtained in a straight-forward 
way and are easy to use, they will hopefully serve a large variety of purposes 
and a wide community in High Energy Physics to obtain interesting new 
results for tree- and loop-level particle decays or resonance cross sections. 
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